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Abstract 

QCD sum rules are used to calculate the q 2 dependence of the nNN coupling 

9ttNn{<1 2 ) in the spacelike region 0.5 GeV 2 < q 2 < 1.5 GeV 2 . We study 

the Borel sum rule for the three point function of one pseudoscalar and two 

nucleon currents up to order four in the operator product expansion. The 

Borel transform is performed with respect to the nucleon momenta, whereas 

the momentum q 2 of the pseudoscalar vertex is kept fixed at spacelike values. 

The results can be well fitted using a monopole form with a cutoff mass of 

about « 800 MeV. 
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I. INTRODUCTION 



The pion nucleoli form factor g n NN{q 2 ) plays a very important, but not less controversial 
role in the framework of ttN and NN dynamics. 

In general, meson baryon form factors are used in one boson exchange potentials (OBEP) 
for the NN force in order to account for the microscopic, i.e. quark and gluon structure of 
the mesons and baryons at a given vertex, and therefore provide a natural cutoff description 
for the interaction potential at short distances. Most of the realistic OBEP fit the data with 
a monopole form for the irNN form factor 



JnNN{q J - 7— 2 , - — T~2 IT U-U 

g-KNNym^) A w - q 2 



with a large monopole mass > 1.3GeV [0]]. On the other hand almost all other hadronic 
theories advocate a much softer form factor with A^ ~ 500 — 950MeV: Chiral soliton and 
quark models for the nucleon the Goldberger-Treiman discrepancy between q 2 = 

and q 2 = m^ 2 ||, general considerations on the structure of the ttNN vertex ||, charge 
exchange reactions |L(J, threshold 7r production [|TT|] , deep inelastic lepton nucleon scattering 
(Sullivan process) p!2]-14|. More recent studies within the OBEP, which either include the 
tc' resonance ||15|| , use a different form of the scalar exchange potential [16|] or consider a 
correlated p7rexchange |T^] also indicate a softer cutoff mass (A^ f& 700 — 800MeV). 

For these reasons it seems highly desirable to perform a calculation of g^NN^ 2 ), which 
does not suffer from the ambiguities of the model and parameterization schemes mentioned 
above and is connected to QCD as closely as possible. 

In this context a quenched lattice QCD calculation has recently been carried out render- 
ing a monopole cutoff mass of A„- = 750MeV Jl8| . 



It is the aim of our work to study this phenomenon within the framework of QCD sum 



rules ||19|| , which have turned out to be a very successful method for calculating hadronic 
properties at intermediate energies (for reviews c.f. pO| , pT| ) without employing the computer 
time consuming lattice gauge calculations. 
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In order to calculate the pion nucleoli coupling constant g n NN = girNNijn 2 ,) itself, i.e. 
for on shell n, one can consider either: (i) the vacuum three point correlator of two nucleon 
and one pseudoscalar meson interpolating fields, which are saturated with resonances in the 



nucleon N and pion it channels on the phenomenological side [EO,EI3-|24| or (ii) the 2 point 



function of 2 nucleon interpolating fields sandwiched between the vacuum and one n state 



and saturating only with N resonances pO|,P3|,P5l . With both methods it was possible to 



obtain rather reasonable results for g n NN, although it should be stated that the uncertainties 
in the Borel analysis are relatively high, especially in method (ii), even if higher order power 
and a s corrections are taken into account |25| . 

If one is interested in the momentum dependence of g^NN^Q 2 ) at intermediate q 2 and 
therefore with off shell 7r, only method (i), which has been successfully used for the calcula- 
tion of meson formfactors p6|-p8|, can be applied. 

Our paper is organized as follows: In section |I| we introduce the 3 point function for 
the NNn vertex and saturate with nucleon intermediate states. In section [TTI| the operator 
product expansion (OPE) is carried out. In section [TV] we perform the Borel analysis. The 
results are discussed and summarized in section [V]. 

II. THE THREE POINT FUNCTION FOR THE vriViV VERTEX 

Our starting point is the 3 point function (Figjl]) 

A(p uP2 ,q) = J d 4 Xl d A x 2 e ip ^e-^ (0^77(2:1) J?(0)^(x 2 )|o) (2.1) 
of a pseudoscalar, charge neutral current 

J 5 °(x) = q(x)i l5 T°q(x) (2.2) 



and two Ioffe nucleon interpolating fields p 



■n{x) = e abc [(u a (x)C lfl u b (x)) l5 Yd c (x)} (2.3) 
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Hereby q = Qj denote the spinors for the quarks with mass mo = m u = ma, a, b, c the color 
indices and C = 27270 the charge conjugation matrix. The three momenta at the vertex are 
related by q = p\ — p 2 . 

Due to restrictions from Lorentz, parity and charge conjugation invariance the 3 point 
function A(pi,p 2 ,q) has the general form 

Mpi,P%,q) = Fi(p\,pl,q 2 )^ + F 2 (pl,pl,q 2 ) #75 + 

F 3 (pl,pl q 2 ) Pl 5 + F 4 (pl,pl q 2 )a, vlb ptp" 2 (2.4) 

where q = pi —P2, P = pi t, P2 . The functions Fi,F 2 ,F i are symmetric and F 3 is antisymmetric 
under the interchange p\ p 2 . 

The matrix element of the pseudoscalar current J5 between on shell nucleon states defines 
the pseudoscalar nucleon formfactor 

(JV(pO| J°(0)|iV(p 2 )) = gp(q 2 )u N (p 2 )z l5 u N ( Pl ) (2.5) 

where ujv(p) denotes a free nucleon spinor in momentum space. 

The pion nucleon coupling constant g n NN is defined by the tiN interaction: 

Ariv = W-kNnNi^ttzN (2.6) 

with the form factor g n NN(q 2 ) For on shell nucleons it is equivalent to use the coupling of 
eq. (|2.6|) or the chirally invariant vector coupling 

By means of the chiral Ward identity 

d^A" = m q(x)i>y 5 Tq(x) (2.8) 
for the axial current A^(x) = q{x)^^^^Tq{x) as well as the PCAC relation 

d^AUx) = mlUir{x) (2.9) 
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one can relate the pseudoscalar nucleon form factor gp(q 2 ) (eq. |2.5| ) to the pion nucleon form 
factor g^NNiq 2 )- 



,rUh = (^)^P, (2.10) 



Therefore saturating the eq.(|2.1|) with nucleon states \N(p)) renders, after continuing to 
Euclidean momenta (q 2 -> -q 2 , p\ -> -p\, p\ -> -p 2 ) 

where A at is the overlap between the Ioffe current r/(x) and a nucleon state 

(0\r){x)\N(p)) = X N e- ipx u N (p) (2.12) 
Using p\ = p\ = Mfq the [. . .] term in eq.( [2.11|) can be cast into the form 



(fa + M N ) l5 {^ 2 + M N ) = (- Pl ■ p 2 + M 2 N ) l5 + ia^ptt&jB + M N ^ 7s (2.13) 
The contribution of the first higher resonance N* with mass to eq.(|2.1|) has the form 



r a i „ - x x (m 2 jA g,NN*{q 2 ) [(fa + M N ) l5 (jk + M* N )] 

[A(p l} p 2 ,q)] N * =i\ N \ N * {——) ^ | _ 2 - 2 ; nj2w 2 ; + JV<->JV 

N 



m / g 2 + m2 (pf + M^) (p| + M 



(2.14) 



III. OPERATOR PRODUCT EXPANSION 



As already stated in eq.( [2.4| ) there are essentially 4 Lorentz structures contained in 
A(pi,p 2 , q). In the following we work with massless quarks, i.e m = 0. From a simple 



dimensional analysis of both sides of the sum rule it is easy to see p0| , |22"| , |2"3"|| that in the sum 
rule for the functions F\ and F4, which both contain an even number of external momenta, 
only operators of even dimension (II, (G 2 ), (qTqqTq), etc.) contribute, whereas the odd 
dimension operators enter with a mass factor tuq and therefore vanish for m — > 0. On the 
other hand for the functions F 2 and F3 containing an odd number of external momenta, 
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only operators of odd dimension ((qq), (qG ■ erg)) contribute. By multiplying both sides of 
the sum rule with jq and taking the trace over the Dirac matrices, we can single out the 
odd dimensional structure, i.e. the sum rule for F2 and F3. Up to order 4 we only have to 
account for the diagrams in Fig.|2| containing the quark condensate (qq) = (uu) = (dd). 
The contribution from the diagrams in Fig.^ja renders: 



[A(pi,P2,q)] 2 a = 16 75 (qq) 4 \ [pIHpD + pI i (pD 

q * 

whereas for the diagrams of Fig.^jb one obtains after a rather lengthy calculation: 



(3.1) 



[A(pi,P2,q)] 



21) 



- 16) 7s (qq) ~ [(pI + pDApIpI q 2 ) + (pi - pl)K + (pl, P 2 2 , q 2 )_ 

- 16) 75 (qq) ^ hi - pl)J(pl,pl Q 2 ) + (pi ~ pl)K„(pl,p 2 , q 2 ) 



( - 16) 7s (qq) 
( - 16) 75 {qq) 



2 

i 

4 

L 

2 



+ 
+ 



\i(p\) + \i{pD + 2% 2 : 
\i{p\) - \i{pD 



+ 



(3.2) 



The expressions I,J,K + and K_ arise from loop integrations and can be most conveniently 
represented as double parameter integrals 



I{P 2 



d A k 



(2tt) 2 (p - kfk 2 



16vr s 



In 



P 

J 2 



(3.3) 



J(p 2 ,pl,q 2 ^ 



and 



d 4 k 



d A k 



(2tt) 2 k 2 (k - Pl ) 2 (k - pa) 



16tt 5 



/o dp L d \ 2 a-x)(i 



o "Jo p 2 (l-A)(l-p)+p|(l-A)p + g 2 Ap(l-p) 



(3.4) 



with 



(2^) 2 P(fc-p 1 ) 2 (A;-p 2 ) 



K ± ( P 2 1 y 2 ,q 2 ) = J o 1 dpJ o 1 d\-^— 



16tt 2 



Ap± (A - 1) 



(3.5) 



(3.6) 



"p? Ap(l - A) + p 2 (l - A)(l - p) + g 2 (l - A)p 
In the UV divergent integral of eq. ( |3.3p we have applied the standard dimensional renormal- 
ization at the renormalization point \x. 
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IV. SUM RULE AND BOREL ANALYSIS 



For the OPE to be valid q 2 has to large, i.e. q 2 ^> Aq CD . In this case we can neglect the 
pion mass in the pole term in eq. (|3.1| ), which is consistent with putting mo = 0. 

In refs. p0| , ^2| , |23|| a sum rule for the pion nucleon coupling g n NN = dnNN^q 2 = 0) has 



been obtained and analyzed by identifying the residua of the pole in eq.( |2.11|) and in the 
OPE contribution from Fig.|2|a (|3.1| ). By introducing 

kg^NNiq 2 ) = g-KNN{q 2 ) - 9ttnn(0) (4.1) 

and analogous expressions for the higher resonance contributions, we therefore can write 
down the sum rule: 



iX » 1 iZt) (p? + Mm + M 2 N ) A9 - NN{q2) + ••• = if^\4A{ P ^q)} 2b 

(4.2) 



Following refs. |2q-p8| we apply the double Borel transform with respect to the nucleon 
momenta p\ and p\ 

B n =B 1 B 2 (4.3) 

2 



r 2 = Ei 

i n t 

and put M\ = Mf = M 2 . The momentum q 2 of the pion is kept fixed at a spacelike value 
q 2 ^> Aq CD , so that both sides of the sum rule depend on M 2 and q 2 now. 

After applying the operator Bu on the OPE side the terms which depend only on either 
p\ or p\ vanish. The Borel transform of the other expressions can be performed by using 
the relations 

\p 2 + sj M 2 

\p 2 + sj M 2 
B{(p 2 ) k e-^ 2 } = {-f±-^ k \a - ±) (4.5) 
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Taking the first resonance with the quantum numbers of the nucleon N* explicitly into 
account in addition to the ground state N we obtain the Borel sum rule 



\2 / " l l 



mlfA M 



N -2-M- 

e m2 



2\m\ 



m J M 6 
'mlfA M N + M* N 



[Ag^ NN (q 2 



+ 



m / 
M 2 ) 2tt 2 



2M 6 



T M* 
N N 

e m 7 e 



q 



where the function g is defined as: 



dp 



x 



+ 2 



1 + p 2 - 4p 



(4.6) 



(4.7) 



;i+p) 3 (i+p) 4 

In order to eliminate the parameters of the higher resonance contribution as far as possi- 
ble we can take the derivative „, \ . on both sides of eq. ( |4.7| ) and substitute the so obtained 
sum rule back into eq.(|4.7|). After employing the Gell-Mann-Oakes-Renner relation 



m nL = -2m (qq) 



(4.8) 



we obtain finally 



with 



B [q 2 ,M 2 ,M 



* 2 
N 



Ag nNN (q 2 



M 2 



M N 4tt 2 X 2 n 



B(q 2 ,M 2 ,M 



* 2 
N 



q 2 \ M* N 2 + M 2 N - 3M 2 



M* N 2 



M 2 N 



+ 



M* N 2 



M 2 N 



M 2 



M 



M 



N 



(4.9) 



Afi 

2 

e 'm t 



(4.10) 



Before performing the Borel analysis, let us comment about the region of q 2 , where our 
method can be assumed to work. If q 2 is too low (q 2 < ^-qcd) the OPE breaks down 
due to higher order power corrections in \. On the other hand we have used the PCAC 
relation ( |2.9[) in order to relate g^NN^q 2 ) to the pseudoscalar nucleon form factor gp(q 2 ) 
(eq. fl2.10p ). This is equivalent to assume pion pole dominance for gp(q 2 ), i.e. saturating 
(N (pi)\J®(0)\N (p 2 )) with a pion state and using the coupling g^NN^q 2 )- This assumption 



likely breaks down if q 2 is of the same magnitude as the first pionic resonance n', i.e. 
q 2 ~ m 2 , ~ 2.2Mjf, because then ir' contributes in the dispersion relation to a similar 
extent than 7r does. 



Keeping this in mind we use for the the Borel analysis q 2 values in the interval 0.5M 



2 < 

N ~ 



q 2 < 1.5M^, where our approach can be assumed to be reliable. Hereby we follow the spirit 
of ref. []30| and treat the resonance mass Mjy as an effective mass, which value is adjusted 



in order to obtain maximal Borel stability. 

From Figs. (fUJi!) we can deduce that for the whole q 2 interval mentioned above an 
effective resonance mass of M^ 2 « 6.0M^ (Mjy « 2.3MeV) gives the largest Borel plateau 



of 0.8M& < M 2 < 1.8M%. 



V. RESULTS AND DISCUSSION 



From eqs.( |4.9] , (4.10| ) we obtain the final expression for the form factor fnNNiq 2 ) defined 
in eq. ( |1 . 1| ) : 



UNN{q 2 ) = 1 - ~TT~ M } 2X 2 B (1 2 > M N 2 ) (5-1) 

M N g^NN^Xjf 



Due to the discussion in the last section, we use M 2 = Mjy and M^ 2 = 6.0M^. For Mjv 
we take the experimental value = 940MeV. In order to avoid uncertainties from the 
sum rule for q^nn we will also use the experimental value q^nn = 13.4 instead of the value 
from the sum rule of refs. p0|j2^j23[| . The parameter is not experimentally known but 



has been determined from various analyses using the nucleon sum rule, the most recent and 



concerning uncertainties obviously most reliable one coming from ref. |31 



Xn = 5.5 • l(T 4 GeV b = 8.07 • 1(T 4 M^ (5.2) 



As we have discussed in the last section formula (|5.2|) cannot be applied for low q 2 . This 
means that we have to interpolate from q 2 « 0.5Mjy to q 2 = 0, where f^NN^) — 1- It is 
interesting to not that the values for f n NN(q 2 ) which one obtains by this procedure at low q 2 
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are practically identical to the ones one would get if one applied eq. ( |5.2|) literally to the low 
q 2 region. This might be purely accidental, but a possible reason could be that higher power 
corrections in^ blowing up the OPE at q 2 — ► are canceled or suppressed after applying 
the double Borel transform B12 ( |4.4| ). In order to answer this question one has to examine 
higher order power corrections to the OPE diagrams in Fig|| which is rather cumbersome 
for a 3 point function consisting of four quark lines and will be postponed to a separate 
analysis. 

Anyhow, from Fig.||we can see that the interpolation of fnNNiq 2 ) from intermediate q 2 , 
where we can apply eq.( |5.2|) , to small q 2 is very smooth. Furthermore in the whole interval 
< q 2 < 2.0M^r the f^NN^q 2 ) of our calculation can be fitted very accurately by a monopole 
form ( [L.l|) with a cutoff mass of A„. = 0.85Mtv = 800MeV. This result is very close to the 
one recently obtained from quenched lattice QCD fl8j , which indicates that both methods 
give a similar description of the ttNN vertex at intermediate q 2 . 
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FIGURES 



FIG. 1. The three point function A(pi,p2,q). 



FIG. 2. Diagrams in the OPE 



FIG. 3. The Borel curve B(M\q\M* N z ) at q* = 0.5M N Z for 3 values of the resonance mass 
M* N 2 . Maximum stability is obtained for M* N 2 = 6.0M N 2 . 



FIG. 4. As in © for q 2 = 1.0M N 2 . 



FIG. 5. As in (Eh for q 2 = 1.5M N 2 



FIG. 6. The fnNN(q 2 ) calculated from eq.(U) (interpolated to small q 2 ) with M 2 = M N 2 
and M^j 2 = 6.0Mjv 2 (full line) compared with a monopole form ( |1 . 1[ ) with a cutoff mass of 
A n = 0.85M N = 800 MeV (dotted line). 
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Monopole Fit 




